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Abstract 
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10~ 5 with a total relative error of 0.8%. This measurement complements the study of S- 
and P-wave hadronic form factors by assigning absolute values to the relative hadronic form 
factors obtained earlier in a simultaneous analysis of the 7T7T scattering lengths conducted on 
the same data sample. The overall form factor normalization f s — 5.705±0.017 oxp ±0.031 cxt 
is obtained with a total relative precision of 0.6%. 
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1 Introduction 



The interest of — > TT + TT~e ± i / decays (denoted K e t± in the following) was recognized many years 
ago at a time when only a handful of such events had been observed pQ. The accumulation of 
a large sample of more than one million of such decays by the NA48/2 experiment has recently 
allowed a very detailed study of the -kit scattering lengths and hadronic form factors [2] . In that 
study, the 1 = and 1 = 2 S-wave scattering lengths have been determined with an improved 
precision comparable to the few percent relative accuracy of the most elaborate theoretical 
predictions [3]. Without the branching ratio value, only relative form factors could be measured, 
giving a full set of values up to a common normalization. 

A new measurement of the K^ A and K~ A decay rates based on the data collected by the 
NA48/2 experiment at the CERN SPS in 2003-2004 is reported here. The event sample is 
about three times larger than the total world sample and has one percent level background 
contamination. A good control of systematic uncertainties, dominated by the external error 
from the normalization mode, allows rate and form factors to be measured with an improved 
precision. These can be used as input to the determination of the Low Energy Constants (LEC) 
of Chiral Perturbation Theory (ChPT) |U [6] and as tests of other theoretical dispersive 
approaches [7]. 

2 The NA48/2 experiment beam and detector 

The NA48/2 experiment, specifically designed for charge asymmetry measurements [8], takes 
advantage of simultaneous K + and K~ beams produced by 400 GeV/c primary CERN SPS 
protons impinging on a 40 cm long beryllium target. Oppositely charged particles, with a central 
momentum of 60 GeV/c and a momentum band of ±3.8% (rms), are selected by two systems 
of dipole magnets with zero total deflection (each system forming an 'achromat'), focusing 
quadrupoles, muon sweepers and collimators. 

At the entrance of the decay volume housed in a 114 m long evacuated vacuum tank, the 
beams contain ~ 3.6 x 10 6 charged kaons per pulse of about 4.5 s duration with a flux ratio 
K + /K~ close to 1.8. Both beams follow the same path in the decay volume: their axes coincide 
within 1 mm, while the transverse size of each beam is about 1 cm. 

The decay volume is followed by a magnetic spectrometer located in a tank filled with 
helium at nearly atmospheric pressure, separated from the vacuum tank by a thin (0.3%Ao) 
Kevlar® window. An aluminum beam pipe of 16 cm outer diameter traversing the centre of the 
spectrometer (and all the following detector elements) allows the undecayed beam particles and 
the muon halo from decays of beam pions to continue their path in vacuum. The spectrometer 
consists of four octagonal drift chambers (DCH), each composed of four staggered double planes 
of sense wires, and located upstream (DCH1-2) and downstream (DCH3-4) of a large aperture 
dipole magnet. The magnet provides a transverse momentum kick Ap = 120 MeV/c to charged 
particles in the horizontal plane. The momentum resolution achieved in the spectrometer is 
a p /p = (1.02 0.044 • p)% {p in GeV/c). 

The spectrometer is followed by a hodoscope (HOD) consisting of two planes of plastic 
scintillator segmented into vertical and horizontal strip-shaped counters (128 in total). The 
HOD surface is logically subdivided into 16 exclusive square regions whose fast signals are used 
to trigger the detector readout on charged track topologies. Its time resolution is ~ 150 ps. 

A liquid krypton electromagnetic calorimeter (LKr), used for particle identification in the 
present analysis, is located behind the HOD. It is an almost homogeneous ionization chamber 
with an active volume of 7 m 3 of liquid krypton, segmented transversally into 13248 projective 
cells, approximately 2x2 cm 2 each, 27Ao deep and with no longitudinal segmentation. The 
energies of electrons and photons are measured with a resolution cje/E = (3.2/y^E © 9.0/E 1 © 
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0.42)% (E in GeV) and the transverse position of isolated showers is measured with a spatial 
resolution a x = a y = (0.42/-/E 0.06) cm. 

The muon veto counter (MUV) is located further downstream. It is composed of three planes 
of plastic scintillator slabs (aligned horizontally in the first and last planes, and vertically in the 
middle plane) read out by photomultipliers at both ends, each preceded by a 0.8 m thick iron 
absorber. The MUV is also preceded by a hadronic calorimeter (not used in this analysis) with 
a total iron thickness of 1.2 m. 

A more detailed description of the NA48 detector and its performances can be found in [9]. 

A dedicated two-level trigger selects and flags the events. At the first level (LI), charged 
track topologies are selected by requiring coincidences of hits in the two HOD planes in at least 
two of the 16 square regions. At the second level (L2), a farm of asynchronous microprocessors 
performs a fast reconstruction of tracks and runs a decision-taking algorithm. This trigger logic 
ensures a very high trigger efficiency for three-track topologies. Inefficiencies are typically a few 
10~ 3 at the first level and a few 10 -2 at the second level (more details can be found in [21 18]). 



3 Branching ratio measurement 

The K e 4 rate is measured relative to the abundant K — > 7r + 7r _ 7r ± normalization channel (de- 
noted Ks n below). As the topologies of the two modes are similar in terms of number of charged 
particles, the two samples are collected concurrently using the same trigger logic and a common 
event selection is considered as far as possible. This leads to partial cancellation of the system- 
atic effects induced by an imperfect kaon beam description, local detector inefficiencies and a 
trigger inefficiency. The K e ^ rate relative to and the K e ^ branching ratio (BR) are obtained 
as: 

rOM/r(K 37r ) = • 4i£» (i) 

and 

BR(K e4 ) = • 4^ • BR(i^), (2) 

where N s , Nb, N n are the numbers of signal, background and normalization candidates (the 
background in the normalization sample is negligible), A s and e s are the geometrical acceptance 
and trigger efficiency for the signal sample, A n and e n are those of the normalization sample. 
The normalization branching ratio value BR(i^3 7r ) = (5.59 ± 0.04)% is the world average as 
computed in [TO] . 

It should be noted that the K~ 4 rate has never been measured. As no difference is expected 
from the rate [JJ, a comparison of the separate measurements of the K + and K~ rates is 
used as a consistency check. 



3.1 Event selection 

The same data sample has been considered in both signal and normalization studies. Given 
that BR(i^3 7r )/BR(X e 4) ~ 1400, a filtering of the data stream and the analysis are performed in 
such a way that the candidates are effectively prescaled by a factor of 100 with a negligible 
error, while the K e t± candidates are not affected, leading to a significant reduction of the data 
volume. 

The analysis of the tttt scattering lengths and form factors presented in [2 J focuses on a sample 
free of hard radiative events at the price of some cuts on the additional photon activity in the 
LKr calorimeter. The present analysis includes radiative events and thus loosens or removes 
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some of the selection cuts which could bias the event counting because of imperfect modeling of 
the photon emission mechanism. 

Common selection 

Three-track vertices (compatible with either K e ± or decay topology), in events satisfying 
the three-track trigger logic conditions, are reconstructed by backward extrapolation of track 
segments from the spectrometer into the decay volume, taking into account the measured stray 
magnetic field in the vacuum tank and multiple scattering. The reconstructed vertex must 
satisfy the following criteria: 

- total charge of the three tracks (called "vertex tracks" in the following) equal to ±1; 

- longitudinal position of the vertex within the fiducial decay volume, 2 to 95 m downstream 
of the final collimator, and its transverse position within 5 cm of the nominal beam axis; 

- vertex tracks consistent in time within 12 ns; no additional in-time tracks present in the 
reconstructed event (see section 13.41 for details) ; 

- all vertex tracks within the DCH, HOD, LKr and MUV geometric acceptances; distance 
between any track and the beam mean position (monitored with decays) in the DCH1 
plane greater than 12 cm for better time-dependent acceptance control; 

- track separations required in the DCH1 and LKr planes (minimum allowed distance 2 cm and 
20 cm respectively) to suppress photon conversions and to ensure efficient particle identification, 
minimizing shower overlaps; 

- distance from the impact point of each vertex track on the LKr plane to the closest inactive 
cell of the calorimeter larger than 2 cm to provide maximum collection of energy deposit; 

— total momentum of the three tracks | YlPi\ below 70 GeV/c; 

- no track-associated signal allowed in at least two planes of the MUV in-time with any vertex 
track (within 10 ns). 

If several vertices satisfy the above conditions, the one with the lowest fit x 2 1S considered. 
Particle identification 

Particle identification criteria are based on the geometric association of an in-time LKr energy 
deposition cluster to a track extrapolated to the calorimeter front face (denoted "associated 
cluster" below). The ratio of energy deposition in the LKr calorimeter to momentum measured 
by the spectrometer (E/p) is used for pion/electron separation. A track is identified as an 
electron (e^) if its momentum is greater than 2.75 GeV/c and it has an associated cluster with 
E/p between 0.9 and 1.1. A track is identified as a pion (vr 1 * 1 ) if its momentum is above 5 GeV/ c 
and it has either no associated cluster or an associated cluster with E/p smaller than 0.8. 

Powerful further suppression of pions mis-identified as electrons within the above conditions 
is obtained by using a discriminant variable which is a linear combination of quantities related 
to shower properties (E/p, radial shower width and energy weighted track-cluster distance at 
LKr front face), and is almost momentum independent. The discriminant variable was trained 
on dedicated track samples to be close to 1 for electron tracks and close to for pion tracks 
faking electron tracks (the discriminant variable performances are shown in section [3"T4"|) . In the 
signal selection, its value is required to be larger than 0.9 for the electron track. 

Signal sample 

The K e 4 candidates are then selected using the following particle identification and kinematic 
criteria: 

- the vertex is composed of a single electron candidate and a pair of oppositely charged pion 
candidates 7r + 7r~; 

- the invariant mass of the three tracks in the 7r + 7r _ 7r :l= hypothesis (M^ n ) and the transverse 
momentum pt relative to the beam axis are outside a half-ellipse centered on the nominal kaon 
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mass [10J and zero pt, with semi-axes of 20 MeV/c 2 and 35 MeV/c, respectively, thus requiring 
a non-zero pt value for the undetected neutrino and rejecting fully reconstructed three-body K^ n 
decays (the K e4 signal loss from this cut is ~ 4.5%, as shown by simulation); 

- the square invariant mass M| in the — > tt^X decay is larger than 0.04 (GeV/c 2 ) 2 to 
reject K — > tt^tt decays with a subsequent 7r° — > e + e~7 decay; 

- the invariant mass of the e + e~ system (assigning an electron mass to the oppositely charged 
pion candidate) is larger than 0.03 GeV/ c 2 to ensure rejection of converted photons and of some 
multi-7r° events (as — > 7r ± 7r 7r ). 

Extra rejection of three-body decays is obtained by reconstructing the kaon momentum under 
the assumption of a four-body decay with an undetected massless neutrino. Imposing energy- 
momentum conservation in the decay and fixing the kaon mass and the beam direction to their 
nominal values, a quadratic equation in the kaon momentum px is obtained. A K e4 candidate is 
accepted if a solution is found in the nominal range between 54 and 66 GeV/ c, allowing a small 
fraction of solutions with negative but close to zero equation discriminant values as observed 
for reconstructed simulated signal events because of non-perfect resolution (in this case, a single 
solution is obtained by setting the equation discriminant to zero). 

A total sample of 1 108 941 K e4 candidates (712 288 K+ and 396 653 K~) were selected 
from a total of ~ 2.5 x 10 10 triggers recorded in 2003-2004. The selection is illustrated in Fig. [1^, 
in the (M^ n ,pt) plane and the reconstruction in Fig. [TJ) by the kaon momentum distribution. 

Pt (MeV/c) 




-60 -40 -20 



20 40 60 80 100 



M 3n -M K (MeV/c 2 



50 52 54 56 58 60 62 64 66 68 70 



PK (GeV/c) 



Figure 1: (a) Reconstructed (M^jPt) plane for the K e4 signal candidates. The elliptic cuts used 
in the signal (solid line) and normalization (dashed line) selections are shown, (b) Reconstructed 
kaon momentum from signal events for data after background subtraction (symbols), simulation 
normalized to data (histogram) and background events (scaled by a factor of 5 to be visible) as 
shaded area. The arrows point to the kaon momentum cuts values. The lower plot is the ratio 
of the two spectra (data/simulation) displayed in the upper plot. 



Background estimate 

The decay is the most significant background source. It contributes either via the decay in 
flight of a single pion (tt^ — > e^v) or mis-identification of a pion as an electron. Only pion decays 
occurring close to the kaon decay vertex or leading to a forward electron and thus consistent 
with a three-track vertex and satisfying the {M^^pt) requirements contribute to the background. 
Other background sources come from — > 7r ± 7r°(7r°) decays with subsequent Dalitz decay of 
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a 7T° (?!£, — )■ e + e~7), an electron mis-identified as a pion, and photon(s) undetected. Such two- 
or three-body decay topologies are very unfavored by the signal selection criteria and contribute 
at sub per mil level. 

Decays violating the AS = AQ rule would appear as "wrong sign electron" (WS) 7r ± 7r ± e =F ^ 
K e 4 candidates and are expected at a neg ligible rate (BR < 1.2 x KT 8 at 90% confidence 
level [10 1). The kinematic distribution of the background events is then to a good approximation 
identical to that of the reconstructed WS candidates multiplied by a factor of 2 as two pions 
from decays can mimic the signal final state while one pion only contributes to the WS 
topology. The uncertainty on this factor of 2 is discussed in section 13.41 

Changing the requirement of a pair of opposite charge pions (7r + 7r _ candidates) in the vertex 
selection to a pair of same charge pions (tt^tt^ candidates) and keeping all other requirements 
unchanged is sufficient to determine the number of events in the WS sample. The distribution 
of the WS K e 4 candidates in the {M^,p t ) plane is displayed in Fig. Efcu Another feature of the 
WS sample is shown in Fig. [2b which displays the reconstructed invariant mass of the dilepton 
system in the signal and WS selections. A peak at the value can be seen as expected from 
decays followed by a pion decay in flight. 

A sample of 5276 K e ^ WS candidates (3276 K + and 2000 K~) has been selected concurrently 
with the signal sample. As K 3n decays are the dominant contributors, the total background 
is then estimated to be 2 x 5276 events, a 0.95% relative contamination to the signal. The 
systematic uncertainty on this quantity is discussed in section 13.41 
p t (MeV/c) 




-60 -40 -20 20 40 60 80 100 0.025 0.05 0.075 0.1 0.125 0.15 0.175 0.2 



M 3w -M K (MeV/c 2 ) M eu (GeV/c 2 ) 

Figure 2: (a) Reconstructed (M3 n ,p t ) plane for the K e/ ± background estimated from WS events. 
The elliptic cuts used in the WS selection (solid line) and normalization (dashed line) selections 
are shown, (b) Reconstructed dilepton invariant mass for K e ^ events. Data are shown as 
symbols, simulation as histogram and background events (scaled by a factor of 5 to be visible) 
as shaded area. 

Normalization sample 

The sample is selected applying the following requirements to events passing the common 
selection: 

— the vertex is required to be composed of three pion ■k + tt^-k ± candidates; 

- the invariant mass of the three tracks in the 7r + 7r _ 7r :l= hypothesis (M 3n ) and the transverse 
momentum pt are inside a half-ellipse (as drawn in Fig. []J,) centered on the kaon mass and zero 
p t , with semi-axes 12 MeV/c 2 and 25 MeV/c, respectively, thus requiring fully reconstructed 
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K^n three-body decays; 

- the total momentum of the three tracks | YlPi\ 1S between 54 and 66 GeV/c. 

The reconstructed invariant mass spectrum is displayed in Fig. [3^. Its measured res- 
olution a^ n = 1.7 MeV/c 2 is in agreement with simulation. The three track momentum sum 
distribution is shown in Fig. [3b. The residual disagreement between data and simulation is 
considered in the systematic uncertainties study. 

The number of prescaled candidates in the signal region is 18.82 x 10 6 (12.09 x 10 6 K + 
and 6.73 x 10 6 K~) with a negligible background. 





60 62 64 

P 3n (GeV/c) 



66 



Figure 3: Distributions of the reconstructed invariant mass (a) and the reconstructed kaon 
momentum (b) for the normalization candidates within the final selection. Data are shown 
as symbols and simulation normalized to data as histograms. The lower plots are the ratios of 
the distributions (data/simulation) displayed in the upper plots. 



3.2 Acceptance calculation 

A detailed GEANT3-based [H] Monte Carlo (MC) simulation is used to compute the accep- 
tances for signal and normalization channels. It includes full detector geometry and material 
description, stray magnetic fields, DCH local inefficiencies and misalignment, LKr local ineffi- 
ciencies, accurate simulation of the kaon beam line, and time variations of the above throughout 
the running period. This simulation is used to achieve a large time-weighted MC production, 
providing a simulated event sample about 20 times larger than the signal sample and 1 /4 of the 
prescaled normalization sample, reproducing the observed flux ratio (K + /K~) ~ 1.8. 

The K e 4 signal channel is generated according to the most precise description of the form 
factors as obtained in [2]. The normalization channel Ks n is well understood in terms of simu- 
lation, being of primary physics interest to NA48/2 [8]. The most precise values of the slopes of 
the Dalitz plot have been implemented |12j . Attraction/repulsion between opposite charge/same 
charge particles and real photon emission using PHOTOS 2.15 |13] are included in both simulations. 

The same selection and reconstruction as described in section 13.11 are applied to the simu- 
lated events except for the trigger and timing cuts. Particle identification cuts related to the 
LKr response are replaced by momentum-dependent efficiencies, obtained from pure samples of 
electron and pion tracks. 

The acceptances averaged over periods with different data taking conditions and over the 
two kaon charges are (18.193 ± 0.004)% and (23.967 ± 0.010)% for the signal and normalization 
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channels, respectively. Due to the detector and beam line being largely charge symmetric by 
design, and due to the data taking conditions, these values are practically identical for K + and 
K~ . The uncertainty on the acceptance calculations due to the limited size of the simulation 
samples (a few 10 -4 relative) is included in the systematic error. 

3.3 Trigger efficiency 

Trigger efficiencies are measured from the data using a minimum bias sample downscaled by 100, 
recorded concurrently with the main analysis data stream. The control trigger condition for the 
first level efficiency requires at least one coincidence of hits in the two planes of the scintillator 
hodoscope (HOD). Control triggers for the second level efficiency consist of first level triggers 
recorded regardless of the second level decision. The overall trigger efficiency is (98.52 ±0.11)% 
in the signal channel and (97.65 ± 0.03)% in the normalization channel. The observed difference 
between the two efficiencies can be explained by the different signal and normalization topologies, 
four-body and three-body decays respectively. Three track events from four-body decays are 
less affected by first level trigger inefficiencies (two tracks in the same HOD square region) and 
by local DCH inefficiencies at the second level trigger [8]. The limited statistics of the available 
control samples have a sizable contribution to the systematic error on the K e $ branching ratio 
measurement. 

3.4 Systematic uncertainties 

A large number of possible effects have been studied and quantified, many of them being upper 
limits. When necessary, a correction is applied to account for any observed bias, and residual 
effects are quoted as systematic uncertainty. The considered contributions are described below. 

Acceptance stability. Many studies have been performed varying in turn the value of 
each cut applied in the common, signal and normalization selections. The maximum deviation 
observed with respect to the value of the reference cut has been quoted as the uncertainty if 
statistically significant. None of the studied contributions are dominant and all are below the 
per mil relative level. Varying the common selection cuts contributes 0.03% to the relative 
systematic uncertainty of BR^-fT^). In the signal selection, the anti-background cuts amount 
to 0.03% while the px cut and residual momentum differences together contribute 0.08%. The 
normalization selection cuts add another 0.08%. Momentum cuts in the particle identification 
contribute 0.05% each when considering electron and pion definitions. 

Time control of the beam geometry and acceptance modeling have been investigated in detail. 
While the acceptances for signal and normalization show a relative variation of ~ 2 percent 
between different data taking periods, related to beam geometry and DCH local inefficiencies, 
their ratio is stable in time. A relative change of 0.07% with respect to the nominal result 
is observed when simulated samples used to compute the acceptances are swapped between 
subsamples of the data before combining them. This value is assigned as an upper limit of 
the systematic uncertainty due to time variation of acceptance and beam geometry. It is fully 
consistent with the variations observed when considering smaller subsamples of the data based 
on kaon beam charge and polarity of the achromat and spectrometer magnets. 

The impact of the limited precision of the measured relative form factors [2] on the signal 
acceptance has also been considered (0.06%). The modeling of the amount of material seen by 
decay particles before the magnet could affect bremsstrahlung emission of additional photons. 
As a result of the absence of explicit cuts on additional LKr activity from photons, the estimated 
4% precision on the simulated material thickness has only a 0.06% impact on the final result. 

All the above uncertainties have been added in quadrature to a total relative contribution 
of 0.18%. 
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Muon vetoing efficiency. The MUV veto requirement in the common selection (Sec- 
tion 13, 1|) is essential in suppressing the background to an acceptable leveQ- Removing the 
MUV requirement in the selection of simulated signal and normalization events increases both 
acceptance values, but their ratio A n /A s remains unchanged within 0.05%, suggesting that the 
rejection of late pion decays in flight does not bias the result. The probability to reconstruct a 
common three-track vertex decreases significantly when one or more pions decay to muons. The 
potential effect of the different number of final state pions in signal and normalization channels 
is therefore minimized at the level of the common selection by requiring the presence of a good 
quality vertex. Stability of the result with respect to muon vetoing is also supported by varying 
the minimum track momentum cut from the nominal 5 GeV/c up to 10 GeV/c, where the 
efficiency of MUV hit reconstruction varies by more than 10% while the relative change of the 
final result is within 0.05%. 

In order to estimate the potential bias from MUV reconstruction algorithm and combinatorial 
effects, the requirement for MUV hit association in the common selection has been modified to 
reject events with hits in all three rather than in at least two MUV planes. The observed 0.16% 
relative difference is conservatively quoted as a systematic uncertainty. 

Accidental activity. Possible accidental activity, either from beam particles or from fake 
tracks (ghost tracks) resulting from DCH hit combinatorial , has been subjected to a dedicated 
study. 

The difference At between the time of each vertex track and their average (the vertex time) 
is required to be within ±6 ns in the common selection. Removing this requirement reveals 
different tails in At distribution for and K e 4 selections suggesting different contributions of 
accidental tracks forming a good vertex. The bias due to this effect is estimated by extrapolation 
from the control regions (-16; -10) ns and (+10; +16) ns to the nominal time window (accounting 
for WS events) as illustrated in Fig. d^,. The subtraction of accidental background leads to a 
—0.12% correction to the result, and the difference between two estimates (based on constant 
and parabolic extrapolation to the central signal window) is quoted as its uncertainty (0.06%). 

An event is rejected if an extra non-ghost track is present within a 6 ns window around 
the three-track vertex time. A conservative estimate of the uncertainty due to the presence of 
accidental tracks not forming the decay vertex is obtained by variation of the above time limit 
up to 35 ns, and it is found to be 0.21%. 

Ghost tracks are not allowed to form an accepted vertex and their presence in addition 
to the considered tracks is ignored in order to avoid bias in K^ w and samples related to 
the different reconstruction probabilities of fake tracks from pions and electrons. The ghost 
track tagging procedure using modified criteria (distance between tracks and quality of track 
reconstruction) has been studied in detail to identify the optimal one, and the residual systematic 
bias is estimated to be 0.04%. 

Particle identification. Different pion identification requirements have been studied, re- 
laxing the E/p condition and recomputing signal and normalization acceptances as described 
in section 13.21 The largest difference to the reference result (0.08%) is quoted as the related 
uncertainty. For the electron identification, varying the cut on the linear discriminant variable 
value between 0.85 and 0.95 (or even removing the cut) changes drastically the background 
contamination (up to a factor of four). Applying the corresponding momentum-dependent effi- 
ciency (Fig. HJd) to the simulation, no bias is observed and a maximum deviation of 0.04% from 
the reference result (obtained for the cut value of 0.90) is observed. The uncertainties from pion 
and electron identification are added quadratically. 

Background estimate. The uncertainty on the scaling factor of 2 used to estimate the 

1 The background increases by a factor of 10 when this requirement is removed. 

2 A ghost track is close to another track by less than 1 cm at the DCH1 plane and is reconstructed with worse 
quality. 
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Figure 4: (a) Normalized distribution of the time difference between the time of each vertex 
track and their average for K e ^ (light shaded area) and Ks n (dark shaded area) candidates. The 
vertical arrows indicate the signal time window. The horizontal arrows indicate the two control 
regions. Only the time difference of the track with the largest At absolute value is plotted, (b) 
Efficiency of the linear discriminant variable as a function of momentum for electron tracks (top 
curves). The solid line corresponds to the cut value of 0.90 applied in this analysis. The lower 
dashed line corresponds to the cut value of 0.95 and the upper dotted line to the cut value of 
0.85. For illustration, the pion track rejection (bottom curve) is displayed for the cut value of 
0.90. 



background based on WS K e ^ candidates has been studied in two ways. In the WS event selec- 
tion, the square invariant mass M\ (Section I3.1j) is computed for both pions and can be used 
to classify further the origin of the event. If the smaller mass squared is above (0.27 GeV/c 2 ) 2 
(corresponding to 2m 7r + with resolution smearing), the event is assigned a factor two weight as 
being K^-like (95.8%), otherwise it is assigned a factor one weight as being K nn o -like (4.2%). 
This rough estimate leads to a factor 1.96. Another estimate, based on a simulated sam- 
ple properly weighted for particle-identification performances, gives a similar ratio Right Sign 
electron/ Wrong Sign electron = (1.94 ± 0.15), confirming the prescription for WS background 
related to pion misidentification. It is also in agreement with the overall factor of (2.0 ±0.3) used 
in the form factor analysis [2]. The ±0.15 uncertainty is propagated to the result as a relative 
uncertainty of 0.07%, based on the above studies. In addition, variation of the background- 
related requirements in the K e ± selection within wide ranges shows excellent stability of the 
result. 

Other sources. Dedicated MC samples simulated without real photon emission were used 
to estimate the impact from radiative effects description. One tenth of the full effect was assigned 
as a modeling uncertainty according to the prescription of |14| . 

Trigger efficiency accuracy is limited by the size of the control samples. The statistical 
precision is quoted as a contribution to the systematic uncertainty. 

Simulated samples used in the acceptance calculations contribute to the systematic un- 
certainty through their statistical precision: this source could be reduced by increasing the 
simulation statistics but already contributes at a very low level. 

Sizable uncertainty arises from the external input BR(i^3 7r ) known experimentally with a 
limited relative precision of 0.72% [TO] . 

Two independent analyses have been compared on a subsample of the data and found to be 
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fully consistent, ensuring the robustness of the procedure. 

The breakdown of the considered systematic uncertainties is displayed in Table [TJ 

Table 1: Summary of the relative corrections applied to the BR(if e 4) value and relative contri- 
butions to the systematic uncertainty. 

Source Correction (%) Contribution (%) 



to BR value to BR uncertainty 
Common to all subsamples 

Acceptance stability 0.18 

Muon vetoing efficiency 0.16 

Accidental activity -0.12 0.21 

Particle identification 0.09 

Background estimate 0.07 

Radiative events modeling 0.08 
Subsample-dependent quoted as a global equivalent 

Trigger efficiency 0.11 

Simulation statistics 0.05 

Total systematics -0.12 0.37 

External error 0.72 



3.5 Results 

The final result is a weighted average of 16 values obtained in eight independent data subsamples 
and both kaon charges. The weight of each input includes error contributions of time-dependent 
statistical origin: event statistics (signal, background and normalization), trigger efficiencies 
and acceptances. This method is more robust against time-dependent conditions than using 
an averaged acceptance and trigger efficiency over the whole data taking period. However, 
due to the careful time-dependent treatment of simulated samples, this potential difference is 
kept below the per mil relative level and is taken into account in the systematic error. Other 
systematic uncertainties (Table [T]) are common to all subsamples and are then quoted as a single 
error on the final result. All input ingredients to the BR(A" e 4) measurement are summarized 
in Table [2j By convention, the uncertainties are assigned to three categories: (i) statistical 
errors from the numbers of K e ^ signal candidates (dominant error), WS data events (used 
for background computation) and normalization events; (ii) subsample-dependent systematic 
uncertainties such as those of trigger efficiencies and acceptance and systematic uncertainties 
common to all subsamples; (iii) the external error related to the uncertainty on the normalization 
mode branching ratio (BR(K 3lT ) = (5.59 ± 0.04)% [ID]). 

The resulting values, including all errors, are found to be: 

r(^ e4 )/r(K 3 ^) = (7.615 ± 0.008 stat ± 0.028 syst ) x 10- 4 (3) 

and 

BR(K e4 ) = (4.257 ± 0.004 stat ± 0.016 syst ± 0.031 ex t) x 10~ 5 , (4) 

where the branching ratio error is dominated by the external uncertainty from the normalization 
mode. The BR(i^" e 4) values obtained for the statistically independent subsamples are shown in 
Fig. El also in perfect agreement with the global value obtained from the whole sample and the 
values measured separately for K + and K~: 

BR(^) = (4.255 ± 0.008) x 10~ 5 , BR(1Q) = (4.261 ± 0.011) x 10~ 5 , 

where the quoted uncertainties include statistical and time-dependent systematic contributions. 
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Table 2: Inputs to the BR(i^ e 4) measurement for K + ,K~ and combined K^. The relative 
contribution of each item to the BR(i^) statistical uncertainty is also shown in the last column. 
Statistical errors on the acceptance and trigger efficiency values (given within parentheses) are 
taken into account in the systematic error (Table H]) and not in the total statistical error given 
in the last row. 





K+ 


K~ 


K± 


BR(K? A ) 
relative error (%) 


Signal events 


712 288 


396 653 


1 108 941 


0.096 


WS events 


3 276 


2 000 


5 276 


0.013 


Normalization events/ 100 


12 090 376 


6 728 544 


18 818 920 


0.023 


K e 4 acceptance (%) 


18.190 


18.197 


18.193 


(0.020) 


acceptance (%) 


23.970 


23.961 


23.967 


(0.041) 


K e 4 trigger efficiency (%) 


98.546 


98.480 


98.523 


(0.108) 


trigger efficiency (%) 


97.634 


97.687 


97.653 


(0.033) 


Total relative statistical error (%) 


0.100 



BR(K e4 ) x 10 5 
4.34 



2003 



2004 



4.3 
4.26 
4.22 
4.18 
4.14 
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Figure 5: K e ^ branching ratio for eight statistically independent samples and both kaon charges. 
The hatched band shows the experimental error (<7 exp = <r s t a t © fsyst)- The total error (shaded 
band) includes the external error. The fit x 2 is 15.85 for 15 degrees of freedom when including 
the time-dependent errors only. Full symbols correspond to K + results and empty symbols to 
K~ results. 
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4 Form factors normalization measurement 



4.1 Formalism 

The K e 4 decay rate (in s _1 ) is described in the five-dimensional space of the Cabibbo- 
Maksymowicz kinematic variables [15) . namely the dipion (S n ) and dilepton (S e ) squared invari- 
ant masses and the three decay angles [9 n (9 e ), the same sign pion (electron) angle in the dipion 
(dilepton) rest frame to the dipion (dilepton) line of flight in the kaon rest frame, and <f>, the 
angle between the dipion and dilepton planes in the kaon rest frame] as: 



dr 5 = nt/ ! -I s 5 p{S-ki S e ) Jh{Sir, S e , cos 9 n , cos e , (/)) dS w dS e cfcos 9 n dcos 9 e d<p, (5) 

where p(S 7T ,S e ) = Xa^ (1 — z e ) is the phase space factor, with X = gA 1 / 2 ^^, S n , S e ), a„ = 
(1 - 4m 2 r /5 7r ) 1 / 2 , z e = m? e /S e and A(o, b, c) = a 2 + b 2 + c 2 - 2(ab + ac + be). The function J 5 , 
using four combinations of F, G, R, H complex hadronic form factors (Fi,i = 1,4), reads |16| : 



J5 = 2(1 — z e ) (Ji + I 2 cos 26 e + ^3 sin 2 9 e ■ cos 2(j) + I4 sin 29 e • cos 4> + h sin 9 e • cos 4> 
+Iq cos # e + I7 sin # e • sin + ig sin 2# e • sin 4> + Ig sin 2 # e • sin 2</>) , 

where 

h = \ ((l + z e )| J Fi| 2 + i(3 + z e )(| J F 2 | 2 + |F 3 | 2 ) sin 2 ^ + 2z e | J F 4 | 2 ) , 

h = ~\ (1 - 2 e ) (l^il 2 " ^(l^l 2 + m 2 ) sin 2 8 W ) , 

J 3 =-i (l-z e )(|F 2 | 2 -|F 3 | 2 ) sin 2 ^, 

h = \ (l-z e )Re(F?F 2 ) sm9 w , 

h = - (Re(F?F 3 )+z e Re(FZF 2 )) sin^, 

h = - (ite(F 2 *F 3 ) sin 2 W - z e Re(F?F 4 )) , 

I 7 =- (Im{FfF 2 ) + z e Im(FlF 3 )) B m8 w , 

/ 8 = ± (l-z e ) Im^Fg) sin^, 

h = - *e) Im(F£F 3 ) sin 2 W . 

In i^ e 4 decays, the electron mass can be neglected (z e = 0) and the terms (1 ± z e ) become 
unity. One should also note that the form factor F^ is always multiplied by z e and thus does not 
contribute to the full expression^]. With this simplification, the complex hadronic form factors 
Fi reduce to: 

F 1 =m 2 K { r iF + aGco&O v ), F 2 = m 2 K ((3G), F 3 = m 2 K (f3 7 H), (6) 

where one uses the three dimensionless complex form factors F, G (axial), H (vector), and three 
dimensionless combinations of the S n and S e invariants: 

a = a n (m 2 K - S n - S e )/2m 2 K , (3 = a^S^Se) 1 / 2 /m 2 K , 7 = X/m 2 K . 

A further partial wave expansion of the form factors F\ 2,3 with respect to the cos 9 n variable 
is considered [16J and, once limited to S- and P-wave terms and assuming the same phase 5 P for 

3 The form factor 7? enters only in the definition of F4 and therefore cannot be addressed in the K e 4 decay 
analysis 
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all P-wave form factors, leads to the expressions for F, G, H: 

F = F s e iSs + F p e^cosfl^ , 



G = G p e iS ? , (7) 
H = Hp e iSp . 

The model-independent analysis in [2] determines simultaneously the four real numbers F s , 
F p , G p , Hp and the phase difference (5 = S s — 5 P ) in bins of S n , S e . 

In presence of electromagnetic interaction, the differential decay rate ([5]) is modified by the 
presence of virtual and real photon emission. This effect is implemented in two steps. First, the 
Coulomb attraction/repulsion between two opposite/same charge particles is considered: 

with ujij = 2naQiQj/ , where a is the fine structure constant, QiQj = —1 for opposite charge 
particles (+1 for same charge particles) and /3y is their relative velocity (in units of c). The 
largest effect comes from the attraction between the two pions at low relative velocity. It depends 
only on the S n variable. The electron (positron) being relativistic, the other attractive/repulsive 
pion-electron terms are constant and their product amounts to 0.9998. Then, the PHOTOS 2.15 
program [13] interfaced to the simulation is used for real photon emission. Its effect is a distortion 
of the kinematic variable distributions and is evaluated on a grid of the 5-dimensional space. 

4.2 Form factor determination 

The K e i decay form factors were extensively studied [2] with the same data sample as used for the 
present analysis and their energy variation described as a series expansion of the dimensionless 
invariants q 2 = (S' 7r /4m^) — 1 and S e /Am 2 . All values have been given relative to a common 
value f s , the S-wave axial vector form factor F s (q 2 = 0, S e = 0) : 



1 + f's/fs q 2 + fslfs q 4 + f'Jfs S e /Aml , 

fp/ fs ) 

g P /fs + g'p/fs q 2 , 

hp/fs ■ (8) 

Integrating dT^ ([5]) over the 5-dimensional space after substituting Fi, F2, F3 ([6|) by their ex- 
pression and measured values (0 [8]) , including radiative effects and leaving out the | V us \ and f s 
constants, the K e 4 branching ratio, inclusive of radiative decays, is obtained as: 

BR(K e4 ) = t k ± • (\V US \ ■ f s ) 2 ■ JdT 5 I (|Ks| • fs) 2 , (9) 

where r^± is the mean lifetime (in seconds). The value of f s is then obtained from the 
measured value of BR(i^T e 4) and the integration result. 

Because of the quadratic dependencies displayed in (|9|), the relative uncertainty on f s is 
only half the relative uncertainty from BR(K e 4), kaon lifetime and kinematic space integral, 
while any relative uncertainty on V us propagates with full size. Contributions are categorized 
as follows: 

— Statistical error stems only from the BR(X e 4) measurement. 
— Systematic uncertainties originate from the BR(X e 4) measurement and the phase space inte- 
gral evaluation. Uncertainties on the integration result when varying each relative form factor 



Fs/fs = 

Fp/fs = 

G p / f s = 

Hp/fs = 
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and energy dependence within ±1<t have been considered. The known large anti-correlations 
between f' s , f'J and g p , g' have been omitted to be conservative. The detailed description of the 
phase shift between S- and P-wave form factors ([7]) has a negligible impact. The robustness of 
the integration method has also been checked against several integration grid definitions. As in 
the branching ratio and relative form factor measurements, one tenth of the full PHOTOS effect 
is assigned as systematic uncertainty on the radiative corrections modeling. 
— External inputs contributing to the f s form factor uncertainty are related to the kaon lifetime 
7~x±, the branching ratio of the normalization decay mode JiR^K^) and the \V US \ value. All 
quantities are taken from [ID] . However, it should be kept in mind that only the product \ V US \ ■ f s 
is accessible by this measurement. 

Table [3] summarizes the error contributions. 

Table 3: Summary of the contributions to the f s form factor uncertainties. 



Source relative contribution (%) 



BR(ET e 4) statistical error 


0.05 


BR(i"Le4) systematic error 


0.20 


Form factor energy dependence (systematic error) 


0.21 


Integration method (systematic error) 


0.02 


Radiative effects in integration (systematic error) 


0.04 


Total experimental error 


0.30 


TSH(K e 4) external error 


0.36 


Kaon lifetime (external error) 


0.08 


\V US \ (external error) 


0.40 


Total error (including external errors) 


0.62 



4.3 Results and discussion 

Given the measured K e ^ branching ratio value and using the world average kaon lifetime 
value (1.2380 ± 0.0021) x 10 -8 s, the measurement of the form factors [2] is now complemented 
by the overall f s normalization: 

\V us \-f s = 1.285 ± 0.001 s tat ± 0.004 syst ± 0.005 ext (10) 
corresponding to f a = 5.705 ± 0.003 sta t ± 0.017 syst ± 0.031 ext (11) 

when using \ V US \ = 0.2252 ± 0.0009 [IP] . 

The obtained f s value and its error can be propagated to all relative form factors now 
displayed with absolute values in Table U] and including an additional normalization error, fully 
correlated over all measured values. 

In addition to the above set of values, it can be of further theoretical interest to quote also 
the S- and P-wave normalized projections of the F\ form factor: 

F 1 / 1 m 2 K = F s e i5 ° + (F p + ah G p ) cos 6 n e iS *, (12) 

namely F s and (F p + a/^ G p ), respectively. As all form factors are obtained in simultaneous fits 
together with the phase difference 5 S — 5 P [2], they exhibit correlations which vary with energy. 
In particular, the fit parameters F p and G p are strongly anti-correlated with a coefficient close 
to unity. The combination G p = G p + 7/a F p shows much less correlation with G p (~ 0.20 at 
most) and is also obtained in the fit. To allow an easy interpretation of the results without the 
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Table 4: Absolute values of the form factor measurements (as defined in dHJ)). There are large 
anti-correlations between f' s , f'J (—0.954) and g p ,g' p (—0.914). The normalization error is fully 
correlated over all form factors. 





= 5.705 ± 


0.003 stat 


± 


0.017 syst 


± 


0.031 cxt 


/. 


= 5.705 ± 


0.035 norm 








f's 


= 0.867 ± 


0.040 stat 


± 


0.029 syst 


± 


0.005 norm 


f'J 


= -0.416 ± 


0.040 sta t 


± 


0.034 syst 


± 


0.003 norm 


f'e 


= 0.388 ± 


0.034 stat 


± 


0.040 syst 


± 


0.002 norm 


fp 


= -0.274 ± 


0.017 sta t 


± 


0.023 sys t 


± 


0.002 norm 


fJp 


= 4.952 ± 


0.057 sta t 


± 


0.057 syst 


± 


0.031 nO rm 


9p 


= 0.508 ± 


0.097 sta t 


± 


0.074 syst 


± 


0.003 norm 


hp 


= -2.271 ± 


0.086 stat 


± 


0.046 syst 


± 


0.014 norm 



explicit description of the fit correlations, the values of F s , G p are given in Table [5] together with 
those of F p , G p and H p . Using f)12|) . the P-wave F\ normalized projection can then be obtained 
as a/7 Gp. It can be noted that for S e = 0, the factor a/7 reduces to a n . 
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Table 5: Absolute values of form factor measurements in ten M vn bins. First error within paren- 
theses is statistical, second is systematic (bin to bin uncorrelated part only). A common relative 
error of 0.62% must be added to each form factor bin by bin measurement, fully correlated over 
all form factor and bin measurements. F s values correspond to the projection of F s (S w ,S e ) on 
the M wn axis. No significant S e dependence has been observed for F p , G p , G p and H p within 
the available statistics. 



Bin 


M n7V barycent' 


sr dimensionless form factors 


number 


(MeV/c 2 ) 


F s 


G p 




1 


286.06 


5.7195(3)(3) 


4.334(74)(19) 




2 


295.95 


5.8123(3)(1) 


4.422(53)(31) 




3 


304.88 


5.8647(3) (2) 


4.550(46)(25) 




4 


313.48 


5.9134(3)(2) 


4.645(41)(23) 




5 


322.02 


5.9496(3)(1) 


4.711(38)(28) 




6 


330.80 


5.9769(3)(1) 


4.767(35)(27) 




7 


340.17 


6.0119(3)(1) 


4.780(34)(30) 




8 


350.94 


6.0354(3)(1) 


4.907(34) (20) 




9 


364.57 


6.0532(3)(1) 


5.019(35)(19) 




10 


389.95 


6.1314(3)(5) 


5.163(36)(21) 


Bin 


dimensionless form factors 


number 




G p 


Hp 


1 




-0.181(67)(15) 


5.053(258)(66) 


-1.795(518)(193) 


2 




-0.324(62)(34) 


5.186(142)(84) 


-2.088(320) ( 77) 


3 




-0.209(60)(33) 


4.941(108)(59) 


-1.995(267)( 98) 


4 




-0.156(58)(32) 


4.896( 91)(51) 


-2.750(246) ( 72) 


5 




-0.366(55)(41) 


5.245( 80)(58) 


-2.045(237) ( 98) 


6 




-0.383(54)(38) 


5.283( 73)(56) 


-2.705(234)( 88) 


7 




-0.218(55)(46) 


5.054( 68)(59) 


-2.203(235)(156) 


8 




-0.302(54)(33) 


5.264( 62)(37) 


-1.856(239)(110) 


9 




-0.309(54)(31) 


5.357( 57)(30) 


-2.096(251)(217) 


10 




-0.264(59)(33) 


5.418( 55)(33) 


-2.865(287)(177) 
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5 Summary 

From a sample of 1.11 x 10 6 K e ^ decay candidates with 0.95% background contamination, the 
branching fraction, inclusive of K e ^ decays, has been measured to be BR(iT e 4) = (4.257 ± 
0.016 exp ± 0.031 ex t) x 10~ 5 using as normalization mode (the experimental error <7 exp is 
the quadratic sum of the statistical <r s tat and systematic <7 syst uncertainties). The relative 0.8% 
precision of the achieved measurement, dominated by the external uncertainty from the normal- 
ization mode, represents a factor of ~ 3 improvement with respect to the world average value, 
BR(if e4 ) = (4.09 ±0.10) x 10~ 5 based on two earlier measurements [17} [T8] . The relative decay 
rate T(K e4 )/T(K 3n ) = (7.615 ± 0.030 cxp ) x 10~ 4 is measured with a 0.4% relative precision, a 
factor of ~ 5 improvement over the current world average value of (7.31 ± 0.16) x 10~ 4 . 

The hadronic form factors that characterize the decay have been evaluated both for absolute 
value and energy dependence. The overall normalization form factor F s (q 2 = 0, S e = 0) has 
been measured with a 0.6% total relative precision as f s = 5.705 ± 0.017 exp ± 0.031 cxt when 
using values of kaon mean lifetime tk and | V us | from [10] , a factor of ~ 2 and 4 improvement 
with respect to the values f s = 5.75 ± 0.08 [H] and f s = 5.59 ± 0.14 [17] obtained by earlier 
experiments. The achieved improved precision on K e ^ rate and form factors brings new inputs to 
further theoretical studies and allows stronger tests of Chiral Perturbation Theory predictions. 

Acknowledgments 

We gratefully acknowledge the CERN SPS accelerator and beam line staff for the excellent 
performance of the beam and the technical staff of the participating institutes for their efforts 
in maintenance and operation of the detector. We enjoyed fruitful discussions about K e 4 form 
factors theoretical descriptions with V. Bernard, G. Colangelo, S. Descotes-Genon, M. Knecht 
and P. Stoffer. 

References 

[1] T.D. Lee and C.S. Wu, Ann. Rev. Nucl. Sci. 16 (1966) 471. 

[2] J.R. Batley et al, Eur. Phys. J. C70 (2010) 635. 

[3] G. Colangelo, J. Gasser, H. Leutwyler, Phys. Lett. B488 (2000) 261. 

[4] G. Amoros, J. Bijnens, P. Talavera, Phys. Lett. B480 (2000) 71. 

[5] G. Amoros, J. Bijnens, P. Talavera, Nucl. Phys. B602 (2001) 87. 

[6] J. Bijnens, I. Jemos, Nucl. Phys. B854 (2012) 631. 

[7] J. Donoghue, J. Gasser, H. Leutwyler, Nucl. Phys. B343 (1990) 341. 

[8] J.R. Batley et a/., Eur. Phys. J. C52 (2007) 875. 

[9] V. Fanti et al, Nucl. Inst. Methods A574 (2007) 433. 

[10] J. Beringer et al. (PDG), Phys. Rev. D86 (2012) 010001. 

[11] GEANT detector description and simulation tool, CERN program library long writeup 
W5013 (1994). 

[12] J.R. Batley et al, Phys. Lett. B649 (2007) 349. 



20 



[13] E. Barberio and Z. Was, PHOTOS, Comp. Phys. Comm. 79 (1994) 291. 

[14] Qingjun Xu and Z. Was, Proceedings of Science, PoS (RADCOR2009) (2009) 071, 
(http:/ / pos.sissa.it / cgi-bin / reader / conf.cgi?confid=92) . 

[15] N. Cabibbo and A. Maksymowicz, Phys. Rev. 137 (1965) B438; Phys. Rev. 168 (1968) 
1926. 

[16] A. Pais and S. Treiman, Phys. Rev. 168 (1968) 1858. 
[17] L. Rosselet et al, Phys. Rev. D15 (1977) 574. 
[18] S. Pislak et al, Phys. Rev. D67 (2003) 072004. 



21 



